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Abstract - We investigate a heat- to charge-current converter consisting of a single-level quantum 
dot coupled to two ferromagnetic metals and one ferromagnetic insulator held at different temper- 
atures. We demonstrate that this nano engine can act as an optimal heat to spin-polarized charge 
current converter in an antiparallel geometry, while it acts as a heat to pure spin current converter 
in the parallel case. We discuss the maximal output power of the device and its efficiency. 



Introduction. — Thermoelectrics in nanostructures 
has generated much interest recently, notably due to the 
strong energy-dependence of transmissions which allows 
large thermoelectric efficiencies [f]. Quantum dots offer 
great versability and unprecedented control over their pa- 
rameters, which has lead to a particular interest in their 
thermoelectric properties [2-7]. Experimentally, thermo- 
electrics has been studied both for open quantum dots [8,9] 
and for Coulomb-blockade dots [10-12]. 

From a theoretical point of view, thermoelectrics in 
multi-terminal setups proves to be interesting as they al- 
low for crossed flows of heat and charge currents. So 
far, two different types of such multi-terminal setups with 
quantum dots have been discussed. In the first, two quan- 
tum dots are capacitively coupled to each other thus al- 
lowing the exchange of energy but not of particles between 
the dots. Both dots are then coupled to hot and cold reser- 
voirs, respectively. In the Coulomb-blockade regime, these 
devices have been shown to be optimal heat to charge con- 
verters that can reach Carnot efficiency [5] . In ref . [6] , the 
scaling with the number of quantum channels of such se- 
tups was obtained for large dots in the form of chaotic 
cavities. 

The second class of devices consists of quantum dots 
coupled to two electronic reservoirs as well as to a phonon 
bath at a different temperature. Both, the case of a single 
dot [13] as well as that of a two-site molecule [14] have 
been discussed in the framework of linear-response ther- 
moelectrics. Additionally, the effect of time-reversal sym- 
metry breaking has been explored for such devices where 
the dot is embedded in an Aharonov-Bohm ring [15]. 

Here, we propose a third class of setups in which a 



single-level quantum dot is coupled to two ferromagnetic 
metals at temperature T& as well as to a third ferro- 
magnetic insulator at temperature Tg, cf. fig. 1. For 
Tb > Te, the quantum dot will absorb magnons from the 
insulator and in turn transfer electrons between the two 
metallic leads. Compared to the coupling of a phonon 
bath to the quantum dot, our setup offers several advan- 
tages. First, magnons provide a more controlled system as 
they are restricted to the ferromagnetic insulator whereas 
phonons are present in any material and, hence, can lead 
to parasitic heat flows. Second, one can drive our quan- 
tum dot heat engine not only by a temperature difference 
between magnons and electrons but also by a nonequi- 
librium magnon distribution that is generated by inject- 
ing magnons into the ferromagnetic insulator electrically 
via the spin Hall effect [16]. Thus, one could operate the 
device at low temperatures where magnons are not ther- 
mally excited and, hence, parasitic phonons are absent as 
well. In addition, the large mean free path of magnons in 
ferromagnetic insulators of several centimeters [17] allows 
to separate the heat engine spatially from the point of 
magnon injection. Third, our setup provides an example 
of a spin caloritronic [18] heat engine that allows to drive 
pure spin currents as well as spin-polarized charge cur- 
rents by thermal gradients. Alternative spintronic heat 
engines have been proposed using nanowires containing 
domain walls [19,20]. In ref. [21] a two-terminal quantum- 
dot setup was discussed that converted a nonequilibrium 
magnon distribution into a directed charge current. 

Finally, while the setup involving Coulomb-coupled 
quantum dots relies on energy-dependent transmissions 
and needs them to be especially tailored to reach the tight- 
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Fig. 1: (Color online) Schematic of the heat-to-current con- 
verter. A quantum dot with a single spin-split level is coupled 
to two ferromagnetic metallic reservoirs at temperature Te and 
an additional ferromagnetic insulator held at a different tem- 
perature Tb- 

coupling limit where Carnot efficiency can be achieved [5], 
our proposal works with energy-independent tunnel cou- 
plings and the tight-coupling limit can be achieved for 
completely polarized metallic electrodes. 

Model. — We consider a single-level quantum dot cou- 
pled to two ferromagnetic metallic reservoirs at tempera- 
ture Te and one ferromagnetic insulator at temperature 
Tb- The Hamiltonian of the system consists of five parts, 

H = 22 Hr + H B + H dot + 22 #tun,r + i?tun,B- (1) 

r r 

The first part describes the two ferromagnetic metals r = 
L,R, 

H r = ^ £rkcrar k(T arkcr- (2) 

Here, a^ ko . creates an electron with momentum k and spin 
it in lead r. Ferromagnetism is incorporated in the sense 
of a simple Stoner model using spin-dependent densities 
of states, p r a{^) = J2k <^( w — £rko-)- In the following, we 
assume these to be energy independent, p ra (io) = p ra - 
The spin-dependence of p ra is characterized by the polar- 
ization p r — (p r -f ~ Pri)/ {Pr^ + Pri) which varies between 
p r = — 1 and p r = 1. The two extreme cases correspond to 
half-metallic ferromagnets with majority spins only while 
p r = describes an unpolarized electrode. 

The ferromagnetic insulator is described in terms of a 
Heisenberg chain 

Hfi = — — Sj • Sj + B z S iz , (3) 
(id) 1 

where indicates summation over nearest neighbors. 

We assume the magnetization of the ferromagnetic ground 



state to point along the +z direction and the weak external 
field B z to be either parallel or antiparallel to the magne- 
tization. We perform a Holstein-Primakoff transformation 
and replace the spin by bosonic operators, 

S i+ - y/2S-blbibi, (4) 

Si- = bly/2S-b\bi, (5) 
S iz = S-b\bi. (6) 

Assuming the length S of the spin to be much larger than 
the average boson number (b\bi), which is a reasonable 
approximation at low temperatures, we can expand the 
square roots in 1/S and retain only the leading order 
terms S i+ = V2Sbi, Si- — V2Sb\. Inserting into eq. (3) 
and performing a Fourier transformation, we obtain the 
Hamiltonian describing the spin wave degrees of freedom, 

H^ = Y J ^b\b CL (7) 
q 

with Lo q = 2 J5 ( 1 — cos q • a) + £? z where a denotes a lattice 
vector of the ferromagnetic insulator. 
The quantum dot is described by 

-Hdot = X! e ^4 c cr + U c\.c^c\ci. (8) 

Here, e a — e =p B/2 denotes the energy of the spin-split 
level. The externally applied magnetic field B is measured 
in units of gp,B where g is the g factor and pb the Bohr 
magneton. The Zeeman splitting determines the energy of 
the magnons that can be absorbed/emitted by the quan- 
tum dot. In the following, we assume B > as only in 
this case, both energy and angular momentum conserva- 
tion can be obeyed when absorbing/emitting a magnon. 
The Coulomb energy U is the energy required to occupy 
the quantum dot with two electrons at the same time. 
In the following discussion, we assume it to be infinitely 
large, U — > oo, such that double occupancy of the dot is 
forbidden. 

The tunneling between the quantum dot and the ferro- 
magnetic metals is given by 

-fftumr = tral^Cv + H.C. (9) 

ktr 

The tunnel matrix elements t r are related to the 
spin-dependent tunnel coupling strengths via T ra = 
27r|t r | 2 p r(T /2. For later convenience, we also introduce the 
total tunnel coupling strength T r = + T,.^. 

Finally, the coupling between the quantum dot and the 
ferromagnetic insulator is given by an exchange interaction 
which in terms of the magnonic operators takes the form 

#tun,B = ^2 M b ti c \ c i + H - c - ( 10 ) 
q 
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Here, we neglected terms of the form (S — &^6 q )cJ.c CT 
arising from the z component of the scalar product be- 
tween the dot and ferromagnetic insulator spin as for 
(&q&q) -C S they can be absorbed into the dot Hamilto- 
nian. We, furthermore, define the spectral weight J(w) = 
27r^ q |j q | 2 (5(w — w q ). In contrast to the tunnel coupling 
strengths, J(uj) will in general have a nontrivial energy- 
dependence that depends, e.g., on the dimension of the 
spin chain [22]. However, as for the following discussion 
this energy-dependence of J{uj) is irrelevant, we will as- 
sume it to be constant and write J(lu) = J. 

Technique. — In order to evaluate the transport prop- 
erties of the quantum dot system, we use a real-time di- 
agrammatic approach [23, 24] in its extensions to systems 
with ferromagnetic [25] and magnonic [26] reservoirs. The 
idea of this approach is to split the system into the strongly 
interacting quantum dot region with only few degrees of 
freedom and the nonintcracting reservoirs with many de- 
grees of freedom. The reservoirs are integrated out using 
Wick's theorem. The remaining quantum dot system is 
then described in terms of its reduced density matrix p. 
In the stationary state, p obeys a master equation of the 
form = W p. The generalized transition rates W can 
be calculated perturbatively in the tunnel coupling to the 
reservoirs. 

As we deal with a setup with collinear magnetizations 
only, the reduced density matrix is diagonal and can be 
written as p — (Po, Pf, Pj,) with the probabilities to find 
the dot being empty (Po) or occupied with spin up (P-f) 
or spin down (Pi). To lowest order in the coupling to the 
reservoirs, the master equation then takes the form 

see eq. (11). 

Here, f^(x) — [exp((x — pi r )/T E ) + denotes the Fermi 
function of lead r with chemical potential p r and f^(x) = 
1 — /+(x). Furthermore, n + (x) = [exp(cc/T B ) — 1] _1 de- 
notes the Bose function and n~(x) = 1 + n + (x). 

In order to calculate the (particle) currents and cur- 
rent correlations, we make use of the framework of full- 
counting statistics adopted for systems that can be de- 
scribed by master equations [27-29]. To this end, we in- 
troduce counting fields x = (XfL, Xjx, XtR> XlR, Xb) for 
electrons with spin up (down) leaving the left (right) 
reservoir and magnons leaving the ferromagnetic insula- 
tor, respectively. We multiply each rate in eq. (11) by 
a factor e lNjsXjs where A^ s is the number of particles of 
type j that left reservoir s — L,R, B in the correspond- 
ing transition to obtain the new matrix W x (For sim- 
plicity, we omit the index j for magnons). The cumu- 
lant generating function &(x) is obtained as the small- 
est eigenvalue of W x . We checked numerically that it 
exhibits the symmetry &(x) — &(~X ~ Ay — Ax), 
where A v = i{eV l4 /T E ,eV L /T E ,eV R /T E ,eV a ./T E ,0) and 
A T = -i((e - B/2)/T E , (e + B/2)/T E , (e - B/2)/T E , (e + 
B/2)/Te,B/Tb), from which fluctuation relations be- 
tween higher-order nonlinear cumulants can be de- 
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Fig. 2: (Color online) Spin-resolved electron and magnon cur- 
rents as a function of the magnon temperature. Parameters 
are J = P, B = 2T E , e = and p = 0.8. 



rived [30]. Interestingly, the above relation holds with- 
out changing the direction of the reservoir magnetiza- 
tions and magnetic fields as such a reversal gives rise 
to a system where all charge (spin) currents flow in the 
same (opposite) direction as for the unreversed system. 
From the cumulant generating function, we obtain the 
currents Ij r and the zero-frequency current correlations 
SiirW = J dt{(i jr (t)i 3 , r ,(0)) - (f jr (t))(J iV (0))] by tak- 
ing derivatives with respect to the counting fields Xjr-> 



In 



(-0 



ae( x ) 



dx. 



Si jr 



H) 2 



ae( x ) 



dXjrXj'r' 



(13) 



(14) 



x=o 



From the spin-resolved electron currents we obtain the 
charge and spin currents as I r — I rcr and 1^ = — 
7 r j., respectively. 

In order to calculate heat currents, we multiply the tran- 
sition rates in the master equation (11) by factors e lEjr ^ jr 
where Ej r is the energy transferred in the corresponding 
tunneling event, measured relative to the chemical poten- 
tial of lead r. As above, we obtain the cumulant generating 
function 6(£) and thereby the heat current Jj r as 



Jj r — ( i) 



dm) 



(15) 



Results. — In the following, we consider a symmetric 
system, i.e., we assume the coupling to the two ferromag- 
netic leads to be of equal strength, Tl = Tr = T. We 
furthermore assume equal polarizations — [prJ = p > 
and a symmetrically applied bias voltage, Vl = — Vr = 
V/2. 

We start by analyzing the spin-resolved electron and 
magnon currents in the absence of a bias voltage. We 
first discuss the basic sequence of transport processes. If 
the dot is initially empty, a spin-up electron enters from 
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one of the two ferromagnetic electrodes. Afterwards, a 
magnon is absorbed, flipping the spin of the electron and 
increasing its energy by B. In a next step, the spin-down 
electron can tunnel out of the dot into the ferromagnetic 
electrodes. Similarly, the reversed process of injecting a 
high-energy spin-down electron, creating a magnon and 
ejecting a spin-up electron is also possible. In thermal 
equilibrium, Tg = Te, both contributions precisely cancel 
and there are no average currents flowing, cf. Fig 2. If 
Tb > Te, the absorption of magnons is more likely than 
their emission. In consequence, we have a magnon cur- 
rent flowing out of the ferromagnetic insulator, a spin-up 
electron current entering the dot and a spin-down current 
leaving the dot. If Tb < Te, the situation is reversed and 
a magnon current flows into the ferromagnetic insulator 
while a spin-down electron current enters the dot and a 
spin-up electron current leaves it, cf. fig. 2. We find this 
picture confirmed by the current cross correlations: They 
are positive for magnon and spin up currents while they 
are negative for magnon and spin down currents. 

In the case of parallel magnetizations, the ratio between 
processes where electrons tunnel in from the left or right 
lead is the same as for processes where they tunnel out 
to the left or right lead. Hence, there is no net charge 
current flowing. Nevertheless, we find a pure spin current 
flowing out of the quantum dot into both leads. As the 
polarization is increased, this current is reduced because 
the tunneling out rate for spin down electrons becomes 
suppressed with 1— p. In the antiparallel case, on which we 
will focus from now on, spin-up electrons enter preferrably 
from the left while spin-down electrons leave preferrably to 
the right. Therefore, we now find a finite, spin-polarized 
charge current through the quantum dot. Analytically, we 
obtain the magnon current 

see eq. (16), 

while the spin-resolved electron currents are given by 
^Lt = ^-Ib, hi = -^tIb, iRt = ~hl and Jrj. = 
—Ijjf. The spin-up current through the left barrier in- 
creases linearly with the polarization p. The spin-down 
current, which has opposite sign, does so as well and van- 
ishes for p = 1. Hence, the larger the polarization, the 
larger the charge current flowing through the system. For 
p = 1 we are in the tight-coupling limit where each ab- 
sorbed magnon gives rise to the transfer of one electron 
from the left to the right electrode. For realistic parame- 
ters of r ~ J ~ 1 GHz, we expect currents of the order of 
a few pA. 

We next turn to the discussion of the Onsager rela- 
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Fig. 3: (Color online) Power and efficiency as a function of the 
applied voltage for different polarizations in the antiparallel 
setup. Parameters are J = T, B = 2T Ey £ = 0, Tb = 2T E . 

tions [31-33] for our setup. Expanding the charge current 
through the left tunnel barrier and the heat current car- 
ried by the magnons to lowest order in the applied bias 
voltage AV and the heat gradient AT = Tg — Te, 

Jl = GAV + LAT, (18) 
Jb = M AV + NAT, (19) 

we find that the Onsager relation M — TeL is fulfilled 
with 

pBT J sech^ 

L= -f- — v . (20) 

4T| (2 J cosh ^ + J cosh %L+T smb. ||J 

Here, we set e = for simplicity. We note, however, the 
Onsager relation is satisfied for any value of e. We remark 
that for the Onsager relation to hold we have to switch 
both the magnetizations and magnetic fields (which for 
our system is equivalent to not switching them at all, cf. 
the discussion above) . Reverting only the direction of the 
magnetic field (B — > —B) leads to an apparent violation 
of the Onsager relation as the coefficients L and M are 
neither even nor odd functions of B. 

We now turn to the power that can be extracted from 
the device in the antiparallel geometry. In order to do 
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+ f ) [f-(e - f ) + n+(B)} - f+(e - f )n+(B) 

f+(e + f ) [r/+(e - f ) - Jn-(B)] - [r + J {l + [2 + f+(e - f )] } n+(B)] 



(16) 



useful work, we need to attach an external load to our 
system. To this end, we apply a finite bias voltage against 
which the thermoelectric current can do work. The output 
power is then simply given by P = JlV. 

In fig. 3, we show the output power as a function of the 
external bias voltage. For V = the output power van- 
ishes as there is no load to work against. For increasing 
bias, the power increases, reaches a maximum and then 
goes to zero again at the stopping potential V s top at which 
7l = 0, i.e., the thermoelectric current is exactly compen- 
sated by the voltage-induced current at this point. We, 
furthermore, see that as the polarization is increased, the 
output power also increases. This is a direct consequence 
of the magnon current being converted more efficiently 
into a charge current for large polarizations which culmi- 
nates in the conversion of one magnon into one transferred 
electron for p — 1. We estimate output powers of ~ 1 fW 
for T - 1 Ghz and B ~ 0.1 meV. 

We next turn to the efficiency of heat to work conversion 
which is defined as the ratio between the output power and 
the input heat. Hence, for Tg > Te where heat flows from 
the insulator into the electronic system, we have rj = P/Jb 
while for Tg > Tb where heat flows from the electronic 
system into the insulator, we have r\ — Pj J2ra Jar- 

The efficiency as a function of the applied bias voltage 
is shown in the lower panel of fig. 3. For p < 1, we see 
that r) = at zero bias and the stopping potential as in 
these cases there is no voltage and no current, respectively, 
and, hence, no output power. In between, the efficiency 
reaches a maximum that increases with the polarization 
p. We note the maximal efficiency in general occurs at a 
different bias voltage than the maximal power. 

The efficiency behaves rather differently for p — 1. Here, 
it again vanishes at V = but then increases linearly with 
the applied bias voltage, reaching Carnot efficiency rjc at 
the stopping potential. This is a consequence of the tight- 
coupling limit reached for p = 1 where heat and charge 
currents are directly proportional to each other [34] . How- 
ever, it is important to realize that for a heat engine at 
Carnot efficiency no work can be extracted. In fact, at 
the stopping potential we have vanishing output power in 
combination with vanishing heat currents. 

The more relevant quantity to judge the efficiency of our 
device is the efficiency at maximum power 77 max p which is 
shown in fig. 4 together with the parameters e, B, V and 
J that maximize the output power for both directions of 
operation. Close to thermal equilibrium, Tb — Te, we 
find that ?7 max p = in agreement with general ther- 

modynamic arguments [35]. We remark that this rela- 
tion usually does not hold for systems whose thermopower 
is not symmetric under time reversal [36]. As we devi- 
ate from Tb = Te, we note that 77 max p becomes larger 
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Fig. 4: (Color online) Efficiency at maximum power as a func- 
tion of Tb (upper panel) and Te (lower panel) for p = 1 in 
the antiparallel geometry. Insets show optimized values of e, 
B and V. The optimization was carried out for the dimension- 
less quantity P/[(J + T)Te/b] for the upper and lower panel, 
respectively. 



than ?7c/2 and reaches Carnot efficiency when the tem- 
perature of the cold bath vanishes. A similar behavior 
was found in Ref. [37] where the nonlinear thermoelec- 
tric efficiency of a quantum dot was discusssed for the 
first time. Interestingly, ?7 max p is larger when the elec- 
tron temperature is larger than the magnon temperature. 
The efficiency at maximum power satisfies the bounds 
< ?7m a xP < ?7c/(2 — rjc)- This is in agreement with 
the findings of Schmiedl and Seifert [38] who gave the up- 
per bound rjss = r/c/(2— rjc) for the efficiency at maximum 
power. 

Conclusions. — We demonstrated that a single-level 
quantum dot coupled to two ferromagnetic metals and 
a ferromagnetic insulator can act as a heat to current 
converter. Depending on the magnetic configuration of 
the metalic reservoirs, either a pure spin current or a 
spin-polarized charge current is generated. We showed 
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the validity of Onsager relations for heat and charge cur- 
rents. We analyzed the maximal output power and demon- 
strated that in the tight-coupling limit the device can 
reach Carnot efficiency Finally, we discussed the effi- 
ciency at maximum power of the magnon-driven quan- 
tum dot heat engine and found that it satisfies the bounds 
r/c/2 < ?7 max p < ?7ss- 
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